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VORTEX MOTION BEHIND A CIRCULAR CYLINDER 

Ludwig Fflppl 
GOttingen 

Presented by S. Fi nsterwalder at the meeting of January 11, 1913 
1. The Equilibrium Condition of the Vortex Pair 

If a circular cylinder in a container of water starts to move from a 
still position, soon after motion commences, two vortices spinning in opposite 
directions form behind the cylinder. The force of the vortices steadily 
increases, since fluid is continually accruing from the boundary layer between 
the edge of the cylinder and the potential current. Meanwhile, the pair of 
vortices moves away from the cylinder, but at a velocity that is small com- 
pared to the velocity through the standing water of the cylinder together 
with the pair of vortices. In the Gottingen Institute for Applied Mechanics, 
the various stages of development of the vortices after their detachment 
were photographed after having been made visible by means of dusting with 
lycopodium powder.*' The photographs reproduced in this article depict the 
vortex pair in an advanced state of development. In the experiment, one 
makes the assumption that the vortex pair which is completely developed is 
almost stationary with respect to the cylinder. 

The observations and experiments cited prompted Herr v. KArmin and 
myself to ask whether there are positions behind a circular cylinder moving 
in still water where the vortices would remain stationary with respect to 
the cylinder. Concurrently, we hoped to deal with the Infinitely large 


1. The photographs were supplied by Herr Rubach who, in his forthcoming 
dissertation, will publish experimental material pertaining to the subject 
that is herein treated theoretically. Herr Rubach has kindly supplied me 
with the accompanying photographs. 

* Numbers in the margin indicate pagination in the foreign text. 
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velocity of the potential current on either end of plates moving through still 

2 

water that likewise generate pairs of vortices. 

I attempted to answer this question by searching for the stationary 
site of the vortex pair behind a circular cylinder of the radius 1. If the 
complex coordinates 

t = £ T *»? 

describe the flow in the plane bound to the cylinder, and the position of 
the vortex is 

^0 = + i% or . *• = *o *%• 

then the complex velocity potential TF can be written with the help of the 
two vortices in the following manner* 


W=-4>+i'P=U(e + i) + iCig } /. 

in which means the velocity in the infinite, and . C the force of the vor- 
tex. 0 is the velocity potential; !F=* a constant aye the streamlines. 

By means of differentiation, one obtains for (l) the complex velocity 


dW 

dt 


h— it’ 



( 2 ) 


In order to derive the velocity of the two vortices which naturally lie /3 
in mirror image across the f-i axis, one substitutes £ with * £„. in (2). 

The term * 



2. This case has already been dealt with years ago by Prof. W. Kutta. An 
extensive characterisation of th« flow can be found In the collection of 
the Mathematical Institute of the Technischen Hochschule in Munich. 
(Footnote by Prof. Finsterwalder.) 
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gives an infinitely large velocity in the center of the individual vortex for 
C = C 0 * which, however, contributes nothing to the displacement velocity 
of the vortex center and so must he floored here. The complex displacement 
velocity of the vortex center is therefore* 



(?) 


If we now leave out the indices 0 as understood and divide the equa- 
tion (3) into real and imaginary parts, then we obtain the two equations 



( 4 ) 


where r- = -f- if . Therefore, the condition for the pair of vortices to 

remain stationary is simply • « = 0 and v = 0, or: 



Since ( appears as a factor on the left side of equation (5b), 
or the v axis is a solution to our problem) that is, a pair of vortices 
whose center points possess the coordinates v and — »; on the »;-? axis can 
remain stationary with respect to the cylinder. The force of the vortex 
belonging to the arbitrarily selected 1 ; can be calculated from ( 5 ) • However, 
this solution for equations (5) has no further interest for us sinoe they 
were not observed during the experiments. Ve seek further solutions for /4 
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the equations (5) or the determinant equation resulting from (5)> 



1)' 5 

1__ 

r* — 1 




( 6 ) 


Since the equation (6) is not changed when the signs of •£ and y are 
changed, the curve must run symmetrically in all four quadrants. Equation (6) 
can be transformed into the following simple form: 

±2 V = (7) 

n 

where the different signs before v belong to the two branches of the curva. 

Equation (7) permits a simple construction of the curve which is de- 
picted in Figure 1 for positive { . It is emphasized that at the point c =1, 
the slope of the curve is less than 4^ and at infinity possesses an 
aysymptote that is angled against the {• axis at less than 30°. The course 
of the curve for negative £ , is naturally the drawn curve mirrored along 
the ?<xis. However, in the Infinite, this branch of the curve has no 
meaning for the direction of the velocity U we have selected r\nd depicted jj> 

in Figure 1 for a certain position of the vortex pair. 


3. The correctness of the transformation of equation (6) to equation (7) 
can be very simply proven by reversing equation (7), which can also be 
written 4 >/*/•* “(«•*— U* , rnd inserting it in equation (6). The left side 
of equation (6) is thus transformed in the following mannert 


1 

1 r* ^ 


Sjl 

r‘ 


1 — 2 1 * 
V-l)» 
1 

>■*- i 


1 c* r 2 1 * 
r* — l" 


2 r®7r* — l) 


(-2 + 2 )'* — )•* + 2 >'* — 1) — (r* — 2 «' fl — 1 + 2 )'»; 
“ 2 1) 


0 . 
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The force of the vortex C, 
determined by the various positions 
of the vortex pair along the curve 
is calculated from (5)» 

( 8 ) 

Accordingly, at a constant 
velocity in the infinite, the 
vortex force C increases as the 
vortex pair lies at increasing 
distances along the curve away from 
the cylinder. 


If we compare the accompanying photographs, then we find by measurement 
that the centers of the vortex pairs lie precisely along the curves deter- 
mined above. From the hydrodynamic events observed when a very long circu- 
lar cylinder is moved from a stationary position, we can make the following 
sketchs A pair of vortices is detached from the rear edge of the cylinder, 
which, with continuously increasing growth of vortex force, approaches our 
curve, and after it has roadbed the curve, continues moving slowly along 
or next to the curve, while the vortex force increases as required by equa- 
tion (8). Now the question is, how does the process develop further. The 
experiment shows that the development of the vortex pair as we have just 
followed it does not continue in the same fashion, but that the develop- 
ment of the flow changes completely. The two vortices abandon their 
symetrical positions with respect to the 1 f- axis. New vortices are formed, 
and now alternately on either side of the cyllnderi a kind of pendulous 
movement is set up in the region behind the cylinder, while the vortices 
detached earlier continue moving at some distance from the cylinder in two 
vortex paths. We now have a flow regime as Investigated by Herr v. 
Karman [\J . The change of the original flow regime is evidently 
caused by an Instability in the current. We shall therefore assign our- 
selves the task of testing the stability of the vortex pair behind the 
circular cylinder. 
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2. Investigation of Stability of the Vortex Pair 

We assume that the flow regime . has developed to the extent that 
the center points of the tiro vortices lie on the curve determined in Section 1; 
we consider, that is, the condition illustrated in the photographs. As was 
already mentioned in Section 1, the vortex pair moves very slowly along the 
curve. Thus, for purposes of investigation of stability, we can view this 
condition as approximating the equilibrium condition. From this equilibrium 
condition, we will displace the vortex pair infinitesimally small distances [ '? 

to see if the vortex pair returns to its equilibrium position. We shall 
conduct the investigation of stability in two partsi First, we will inves- 
tigate stability when displacement from the equilibrium position has a mirror 
image symmetry} then, we will investigate asymmetric displacement. 

To begin with the first case, we shall imagine the vortices displaced 
from their equilibrium position infinitesimally small amounts + a parallel to 
the f- axis and -j -ft or—/? parallel to the »?- axis. Since the two vortices 
thus remain parallel to the , 4 - axis, we can use the equations (4a) and (4b) 
derived for this condition, for the velocities u '=* parallel to the 
axis and a — ~ parallel to the »/- axis. If in equations (4) one sub- 

k u r ^ 5 

stitutes v* with ; + « and >) with »/+/? and develops the left side of equations 
(4) with the small quantities a and ft , retaining only the first potentials 
of a and ft , one obtains, taking the equations (5a) and (5b) into conside- 
ration (which cause the finite terms which appear during development to 
disappear by themselves) we have the following stability equations: 


where 


(a) 

(b) 




A 


X 



(4*;* +1); 


8 ,**»/ 


U 

r*(r* — 1)’ 


3 

Y 



(9) 


( 10 ) 
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If one eliminates one of the two dependent « or ft from the two first 
order linear simultaneous differential equations (9a) and (9b), one obtains /8 
the following second order differential equation which is valid for. ft as 
well as for a , i 


*£ + (r- A) + {BX — A Y)a = 0, 


(U) 


whose general solution is 


where 


a =* C x e'*' -f C t 


/•\ t 


y~a 


± l V{7 — A)- ~ 4 (BX - Tt). 


( 12 ) 


(13) 


The conditions for stability are thus the two inequalities! 

(a) Y> A 

(b) BX~AY>0, (1*0 

which are met, as one can prove by substituting A, 23, X, Y with values from 

( 10 ). 


The stability of the vortex pair for mirror image symmetrical displace- 
ment from the equilibrium position is thus proven. If one imagines a wall 
in place of the ! {- axis, so that fluid is found only on one side, then the 
single vortex is completely stable with respect to infinitesimal displace- 
ment from the equilibrium position, since the other three vortices are 
created as mirror images on the wall or on the cylinder as the case may be, 
and the arbitrary displacement must remain mirror-image symmetrical. Experi- 
ments with separating walls behind the cylinder have not as yet been dealt 
with. 


Investigation of stability of the vortex pair with regard to asaymetrical 
displacement from the equilibrium position remains to be discussed. If the 
coordinates of the equilibrium position of the first vortex are . , and for 


7 
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the second vortex £' — v> then the corresponding coordinates for the dis- 
placed vortex pair are;+«, >] + P or f — a, — »? + /? where a and ft axe again /9 

taken to mean infinitesimally small distances. In this case, we cannot use 
the equations (4) as above, since, like equations (l) and (3)» they are valid 
only for mirror image symmetrical positions of the vortices. However, if in 
general, we characterize the coordinates oi the vortex centers with C, =£, +*» 7 , . 
and = £, -f- then by similar considerations that yield equations (l) 
and (3), we easily find for the complex displacement velocity of the first 
vortexi 



where again £ represents the conjugated imaginary value of f . 


If we separate equation (15) into real and imaginary parts, substitute £, 
with | + a, 1?1 With Tj + P, i s with £ — a and with — q + P and 

solve again for a and p, then we obtain as above the two stability equations* 


(.) A-a + 

. (b) 


( 16 ) 



( 1 ?) 


Just as we derived the equation (11) from the two differential equa- 
tions (9)» *• can use the two differential equations (16) to derive the 
following simplified second order differential equation as a consequence of 
A'-Y' t 


yp+tvr-A ' •).-!), 


( 18 ) 


whose general solution is 


M 


8 



© 




,v 


a ■■ 
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with 
so that 


A'*— B‘X'<0 
A,. 2 «±1/^ ,S ~J5'X', * 


(19) 

( 20 ) 


.is the condition of stability. By inserting the values of * A', B‘, X ‘ 
from (17) • one can easily prove that the inequality (20) is not met for any 
point on the curve, so that the unstabilityof our vortex pair is proved. The two 
root values 


* A, = VA 1 *- B‘X < and A, = _ V A* — B'X‘ * 

belong to the two principal oscillations of the vortex pair for oiymmetricol 
displacements. The value corresponding to \ gives rise to instability, while 
the principal oscillation corresponding to A, represents a dampened move- 
ment. From the first principal oscillation, one obtains from (16) the rela- 
tionship 

A'-VA‘ i -B'X l 


(») 


t 

a 


B' 


and for the latter, the stable motion, 

fi A‘+]/A !r ^~B r X l 


( 21 ) 


(b) 


B 1 


However, siqce according to (17) L"X' is always negative, B‘ always 
positive, and X' 'always negative, one sees that in the case of unstable dis- 
placement " is positive, and in the stable case' jj is negative, and taken 
absolutely, is greater than the first case. If one uses the above formulas 

if 

(21) for ' oi iht specific position of the vortex pair in the photographs, /ii 
one obtains, after simple calculation* 




2 f or unstable displacement 


i 

a 


3.7 for stable displacement. 


The change of the observed current formation which was disoussfc at the 
end of Section 1 thus has Its more fundamental source in the instability >f one 
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of the four principal oscillations of the vortex pair. Indeed, one observes 
that soon after the vortex rair has reached our curve and starts to move very 
a -wly, there is a noticeable change in the hydrodynamic process which is 
cais d by the instability just calculated. However, this instability does not 
appear as a sudden asymmetric displacement of the potential of the two 
vortices as one might expect from the preceding calculation} rather, before 
the its table displacement commences, both vortices lose their characteristics 
as potential vortices, where a spinning nucleus forms in the center of the 
vortex which is deformed into an elipse and occasionally also divides. Only 
after this deformation of the vortex does the displacement referred to above 
occur, which starts the formation of the two parallel vortex streams as 
described by KArmAn. 


The Resistance of the Cylinder 


We shall propose the question! Is the resistance of the moving cylinder 
theoretically explained by the observed movement of the vortex pair away 
from the cylinder? For this purpose, we shall imagine the cylinder moving 
with uniform motion and shall calculate the chronological change of the total 


impulse caused by the displacement of the vortex pair with respect to the 
cylinder. If the vortex pair remained stationary with respect to the uni- 
formly moving cylinder, then the resistance of the cylinder would be zero. 
Indeed, however, one observes that the vortices grow and concurrently move 
on. We shall mathematically grasp this process by assuming that there is a 
potential flow outside the vortex nuclei, but that we nevertheless think of 
the vortex as chronologically growing. The vortex pair is at this time 
moving along the curve determined in Section 1. The total impulse is divi- 
ded into two partsi the impulse of the potential flow around +he circular 
cylinder and that of the vortex pair. The first part is U- rt- o, when Q 
Is the density of the fluid, \U the velocity of the cylinder, and the circu- 


lar radius is 1 [zj , The second part stems from the two solids and the two 
mirrored vortices. The complex potential of these four vortices with the 


coordinates I i, (see Figure 2) 1st 

TV- $ + »!?-. Ci lg 




( 22 ) 
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so that 


(a) 
' (b) 


4> =b C(a — /5) 


( 23 ) 



Figure 2 


The X components of the impulse 
of the total fluid are 

F^g^udxdy^Q^^dxdy' (24) 

where the integration extends over that 

part of the plane which is filled with 
4 

fluid. The integration in terms of 
can he carried out in the last double 
integral, and taking (23a) into con- 
sideration, one obtains 


— Cgfia — f})dij. * (25) 
(circle) 


Zl3 


a here has the meaning indicated in Figure 2, and the integral extends 
over the circumference of the unit circle. The latter Integral, which can 
also be written by substitution of y — sin y , thusi 

2.T 

J «« £(a — /?) cos y dy 

can most simply be solved by means of complex Integration. Ve separate J 
into its two parts which stem from its two pairs of vorticesi 

«'i a cos ydy, J,** ^/tcosydy, 

so that 

and to determine we calculate the complex integral 


4. Correspondingly, one obtains for the Y components the impulse Q-t/ffed* dy, 
However, one can immediately see from reasons of symmetry that V-o . 


11 







ORIGINAL PrtGw s- 
OF POOR QUALITY 


* - f *$=&; * * * - h%=£* <-«*» 


.1 


so that 


S.t 


J\ + J\ — 2 Jig ^~~rr r cos r*r- 


The mearung of the quantities <F and r can be seen in Figure 2. The 
imaginary part is doubled. The two integrals J\ and J\ can, 
however, be easily calculated if one introduces C — ? r as an integration 
variable into the first integral and £ into the second; then: 


( circle) ( circle) 

Since and re~ u r lie outside the path of the integration, then 

•J\ = o . in co^tr^st, J] is different from zero in that 

Jigd — 

(circle) 


2 1.7 


so that 


(circle) 


2 n , . . . 4 jti . 

J\ = — <*r) — — sm <p. 


Thus is 


2 71 

</, = — sin <p. 


The same method is used to calculate J t , and the same formulas are valid 
if one simply substitutes r with - • Correspondingly, one obtains 

J t » 2.7r sin </> , 


so that 


or 


<7 = J, — <7 S = — 2 .7 sin ^ 


«7 ■= — 2 .7 a 


since 


sin ■*«, 

as follows from Figure 2. Equation (25) » which gives the part stemming 
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from the vortices is now 


P = ACaxo -f- 2C«xq = 6C?«.to. 


( 26 ) 


If one now includes the impulse of the potential flow, then one obtains the 
total impulses 


G = Utiq + 6 Cang. 

( 2 ?) 

The fluid pressure on the cylinder is determined by a momentary change of * C ; 
and at constant velocity is 


dG . d(Ca) 


( 28 ) 


Since according to equation (8), can be expressed by the momentary posi- 
tion of the vortex on our curve, then according to (28), the pressure ex- 
perienced by the cylinder at constant velocity can be derived from the 
position and the momentary change of position of the vortex pair. 


Great difficulties stand in the way of calculating the resistance on 
the basis of the observed movement of the vortex pair along our curve. This 
is largely because of the insfabllltyof the vortex proven in Section 2, where 
soon after the vortex pair has reached the curve, it leaves it again as a 
consequence of an unstable disturbance. If, however, two vortex streams have 
formed behind the cylinder as per the vortex formation of v. KArmdn. then, 
as v. KArradn ffj has shown, the resistance of the body can be determined with 
satisfying accuracy from observation of the displacement velocity of the vor- 
tex, as well as the distance between the vortex streams. 


4. The Flow Around the Infinitely Long , Flat Plate 

The observations made in Sections 1 through 3 for a circular cylinder 
suggest expanding the concepts to the motion of a flat plate. We will assume 
that the direction of movement is perpendicular to the plate. It is very 
easy to prove that behind a moving plate there must also be a geometric loca- 
tion for the position of the twovortices turning against the plate with 
opposite rotation. The reason is that if we form the previously considered 


13 



ORIGINAL PAGE IS 
OF POOR QUALITY 


flow around the unit circle of the * plane together with the vortex pair 
through 



( 29 ) 


conforming to e plane of the plate of length 2 and form the derivation of 
the complex velocity potential for the flow in the e plane, which is ob- 
tained from equation (l) by means of elimination of 4 using equation (29), 
then the condition for the vortex pair to remain stationary behind the plate 
is 


dir dvr d: __ 

ds dt'dz~ 0 ' 


(30) 


However, since everywhere in the layer of the z plane -that the physical 

(1 * 

condition jz-r° occurs, equation (30) simplifies to 


dTF 

H=°' 

However, this condition is the same as those expressed in equations (5) 
of Section 1. That is, we obtain the desired geometric location when we ex- 
press the curve determined through equation (7) in the coordinates s = x + iy 
with the help of equation (29) . 

I will not carry out this transformation; rather, I shall deal with the 
most interesting matter of whether there exist positions of the vortex pair [y? 
behind the plate that effect finite velocities at either end of the plate. 

The calculation shows that this condition is met only for the position of 
the vortex pair in the infinite, which is naturally out of the question for 
the explanation of the flow formation. Observations of the flow accompanying 
a moving plate show that similar to a cylinder, soon after movement commences, 
a vortex pair is detached from the rear edge of the plate, which then steadily 
grovs and moves away from the plate. 

If one uses instead of a plate, a bowl with a cross section of a circu- 
lar arc and moves it with the convex side forward, then my theory, as 
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Professor Prandtl pointed out to me, predicts finite velocities at either end 
of the bowl, for a certain position of the vortex pair in the finite and de- 
pending on the flexure of the bowl. I will not elaborate any further on this 
question as no photographs of this flow configuration are at hand. 



Figure 3 



Figure 4 
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